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I Orbital angular momentum (Jo), systemic mass (M) and orbital period 

r~| . (P) distributions of chroniospherically active binaries (CAB) and W Ursae 

^P^i Majoris (W UMa) systems were investigated. The diagrams of log Jo — log P, 

O 

^ ■ log M - log P and log Jo - log M were formed from 119 CAB and 102 W UMa 

^ I stars. The log Jo — logil/ diagram is found to be most meaningful in demon- 

strating dynamical evolution of binary star orbits. A slightly curved borderline 

•i-H . 

. (contact border) separating the detached and the contact systems was discov- 

$^ ' 

' ered on the log Jo — log M diagram. Since orbital size (a) and period (P) of 

binaries are determined by their current Jo, M and mass ratio q, the rates of 
orbital angular momentum loss {d\ogJo/dt) and mass loss {d\ogM/dt) are 
primary parameters to determine the direction and the speed of the dynamical 
evolution. A detached system becomes a contact system if its own dynamical 
evolution enables it to pass the contact border on the log Jo — log M diagram. 
Evolution of q for a mass loosing detached system is unknown unless mass 
loss rate for each component is known. Assuming q is constant in the first 
approximation and using the mean decreasing rates of Jo and M from the 
kinematical ages of CAB stars, it has been predicted that 11, 23 and 39 cent 
of current CAB stars would transform to W UMa systems if their nuclear 
evolution permits them to live 2, 4 and 6 Gyrs respectively. 
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1 INTRODUCTION 



The angular momentum loss in close binaries is known to be replenished from the reservoir 
of orbital angular momentum (QAM) of the system by tidal locking. Since tidal locking 
operates as a mechanism to draw angular momentum from the orbit, the orbits of spin-orbit 
coupled binaries are then forced to shrink although actual angular momentum loss occurs 



at one or bo 



^h components due to magnetica lly driven winds, also called magnetic braking 



( Schatzman 



Kraft 



to be decreasing. Authors s uch. 



Vilhu fc Rahunen I (11 98(1) 



Maceroni fc van't Veer 



Meste] 



(1991 



nski 



19681) . Shrinking orbits, t 



l en, require orbital period s 



Okamo 



Meste] 



19961 ): IStepien 



(199a); 



va.n't Veer 



Guinan Rradstreet 



Demircan 



mi)] 



(19991) all believed this 



mechanism is a main route to form W UMa systems from the binaries initially detached 
with comparable periods. Showing high level of chromospheric and magnetic activity, RS 
CVn like objects including binary BY Dra type systems, shortly called CAB, are primary 
candidates to be the progenitors of the contact systems. 



Karatas et al. I (|2004 ) recently presented observational evidences of mass loss and orbital 



period decrease by comparing the total mass and period histograms between the kinemati- 
cally young (age 0.95 Gyr) and older (age 3.86 Gyr) CAB sub samples. Secular orbital period 
decrease among CAB is indicated further by the fact that the shor ter orbital period svs - 
tems are older than the longer orbital period systems (see Table 5 of iKaratas et al. 1120041) 



Bilir et al. 



Moreo ver, 5.47 Gyr kinematical age for the field W UMa stars determined by 
( 20051 ) appears consistent with the scenario. Consequently, the mean difference of 1.61-Gyr 
age between field W UMa and field CAB could be interpreted as a mean lifetime of contact 
stages. Not all con tact systems have to be formed from detached progenitors. This is because 



Bilir et al. 



( 20051 ) have found that about 20 per cent of contact binary sample have ages 
less than 0.6 Gyr. Such a young age, being much less than the mean life time (1.61-Gyr) of 
contact stages, does not permit pre-contact detached phases. A Pre-contact detached stage 
usually has a duration up to several Gyr, e.g. 3.86 Gyr is an average. Therefore, those young 
W UMa must have been formed directly in the beginning of the main sequence or during 
the pre-main sequence contraction phase. Apparently, both mechanisms with and without 
pre-contact phases actively operate. 



Recently, 



Demircan et al. 



( 20061 ) have determined a first order relative decreasing rates 
of 0AM, systemic mass and orbital period for the detached CAB as Jo/ Jo = —3.48 x 
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10-10 yr-\ 



M/M 



-1.30 X 10~^° yr-^ and P/P = -3.96 x 10"^° yr'^ by using kine- 



matical ages of sub samples formed according to 0AM ranges. Mass loss, 0AM loss and 
orbital period decrease driving one another appear to be confirmed by observational data. 
Nevertheless, it has been known that an isotropic mass loss from the surface(s) of one/or 
both component(s), despite it mean s 0AM loss, forces the orbi t to be bigger and make s the 



orbital period larger (jPringle 



19851 ). It has also been shown by lDemircan et al. I ()2nnfir ) that 



not all rates of 0AM loss would cause an orbit to shrink. Only if there is an amplification 
mechanism, such as magnetic loops forcing magnetized plasma to co-rotate, so mass loss 
at the Alfen radius carries away more angular momentum. If the amplification parameter, 
A = {dJo/dM)/{Jo/M), is bigger than 5/3, then the binary orbit would s hrink and the or 
bital p eriod would decrease. With the decreasing parameters given above, iDemircan et al. 
( 20061 ) found A = 2.68 as a first order observational mean. 

Although, shrinking orbits and decreasing orbital periods are confirmed observationally 
and mean decreasing rates are available to support dynamical evolution, the details of the 
actual process are still unknown. Therefore, in this study, we investigate orbital angular 
momentum (Jo), systemic mass (M) and orbital period (P) distributions of CAB and W 
UMa systems together. 



2 DATA 



Table 1 contains 119 CAB stars from 



Karatas et al. 



( 20041 ) with basic data (orbital and 



physical). Columns are organized as order number, name, spectral type, stage of evolution, 
total mass (M = M1 + M2), mass ratio {q = M2/M1 < 1), orbital period (P), primary radius 
(Pi), secondary radius (P2), logarithm of 0AM (Jo) and spin AM (Jspin) as a fraction of 
0AM. Rows are organized as first 53 stars are possible members of young moving groups 
(MG) and the rest 66 stars are field systems. Furthermore, rows are sorted in increasing 
orbital periods. 

CAB stars in Table 1 had been classified into three stages o f evolution G, SG and MS 



Karatas et al. 



( 20041 ) for 



(column 4) according to their luminosity classes. One can see 
details. Simply, G group contains systems at least one component being a giant. SG group 
contains systems with at least one subgiant but no giants. MS marks the systems with both 
components are on the main sequence. Systems with unknown secondaries were classified 
according to primaries. 
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_A mean kinematical age of 3.86 Gyr was assigned to field CAB stars by 



Karatas et al. 



( 20041 ) from their galactic space velocity dispersions. The field CAB group contains younger 
and older stars as a mixture. But, the o ther CAB group (MG) contains only the young 
ones. One may see iKaratas et al. I (j200J) for the details how MG systems were selected 
from the common CAB. The ages of MG groups are known as open cluster ages by the 
turn off point from the main-sequence. Consequently, a pre-determined age of a moving 
group can be assigned to binaries which were found to be possible members according to 
their space velocity vectors. Therefore, unlike the field stars with various ages, the MG stars 



are homogene ous with a sing. 



considered by 



Karatas et al. 



Gyr age (see Table 3 of 



e age c orresponding to each MG. Among the five MG groups 



( 2004f). t 



le Hy ades super cluster is the oldest one with 0.6 



Karatas et al. 



200J). Therefore, MG syste ms are consider ed to be 



1995 


) and summarized by 
^able 2 contains 102 ^A 


Montes et al. 


( 


2(1 


Ola 


b 




1 


/ UMa stars from 


Bilir et al. 



Eggen 



198S 



( 20051 ) with basic data (orbital and 



physical). Columns are organized the same as Table 1 with only one exception that column 4 
contains the type rather than the evolutionary stage. There are tw o types of W UMa systems 



according to the light and the velocity curves (jBinnendiik 



19701 ). W-type systems are those 



such as W UMa itself, in which the hotter component (the star eclipsed at the primary 
minimum) is smaller and less massive. Rows are organized as first 26 with MG designations 
and the rest 76 stars are field systems. As in Table 1, rows are sorted in increasing orbital 
periods. 

The spin AM (Jspm) of binaries is less well known than the orbital AM because the radii of 
gyration are uncertain. They are needed if one wants to study total AM. They are computed 
for t his study just to compare them to 0AM. An analytical approximation and coefficient 
from IClaret fc Gimenez I ()l990l ) for main-sequence stars were used in computing assuming 
synchronous rotations whenever the radii of both components are available. However, since 
their approximation is for the mass range between 0.60 and 25 Mq, for some CAB and W 
UMa stars with component masses less than 0.60 Mq, Jspm has not been computed. On 
the other hand, the computed Jspin for evolved stars has to be considered as upper hmits 
because the moment of inertia is known to be decreasing by evolution. Those rough values 
of the total spin AM are displayed as fractions of 0AM (as Jspin/ Jo) in Table 1 and Table 
2. It has been found that total spin AM of all systems (CAB and W UMa) usually less than 
one tenth of Jo- There are only nine systems, which are all contact systems (AW UMa, e 
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Cra, FG Hya, OU Ser, TZ Boo, TV Mus, V410 Aur, V776 Gas, XY Boo) with total Spin 
AM higher than one tenth but less a quarter of Jo- In average, W UMa systems have nearly 
10 times more spin AM than those of CAB stars. 

In this study, we consider orbital dynamics, e.g. 0AM changes. If AM transfer occurs 
between the orbit and the spinning components due to tidal interactions, this study considers 
it as an instability in the orbit. Any AM transfered from the orbit to the component stars 
will be sensed by the orbit as an QAM loss and vice versa. For the rest of this study, Jq 
strictly mean QAM, one must not confuse it with the total AM which includes both orbital 
and spin AM existing in a binary. 



3 ORBITAL DYNAMICS 

The most basic definition of 0AM (Jo) can be given as 

where / — (j^q^)a^ = (jj^^^Ma"^ is moment of inertia and = 27r/P is angular speed 
for an orbital motion, thus Jg — I^- 

Jo and M are two basic physical quantities which determine a unique period (P) and a 
unique size for the orbit as 

(l + g)6 27r J! n + qlJl 

where the mass ratio (g = M2/M1 < 1) can be considered as an auxiliary parameter used in 
the definition of Jo- The size a = ai + a2 represents the semi-major axis of a relative orbit 
of one star around the other. Stability of an orbit {dP — 0, da — 0) requires Jo and M to 
be constant {dJo — 0, dM — 0) provided with no mass transfer {dq — 0). If there is no mass 
transfer, which must be true for detached binaries, it is obvious that 0AM loss will cause 
an orbit to reduce its period and size. On the contrary, mass loss has an affect of increasing 
the period and the size. Logarithmic derivatives of (2) give all possible relative changes as 



dP 1 — qdq dJo dM da r,^ — (ldq dJo dM 

Because M has higher power than Jo, the affect of mass loss would dominate. For example. 
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in the case of same relative changes of 0AM and mass (isotropic stellar winds, if dq = 
dJo/Jo = dM/M according to (1), then 



0), 



dP _ dJo _ dM da _ dJo _ dM 

which means mass loss and corresponding 0AM loss will have a net effect on the orbit to 
increase both the period and the size. However, there could be additional causes to increase 
relative 0AM loss, e.g. 0AM loss of gravity waves, or stellar encounters in the galactic 
space, or a third body in an eccentric orbit around the binary system, or existence of an 
amplification mechanism such as in some tidally locked binaries, where tidal interactions 
transfer 0AM to spinning components and transferred AM is lost at the Alfven radius. 
After considering all possibilities, one has to compare the grant total relative 0AM loss to 
the relative mass loss which could be expressed by a parameter 6 defined as 

(5) 

which can be called dynamical parameter because dynamical respond of the orbit depends 
on the value of 6. It should be noted that 6 is equal to the amplification parameter A = 



{dJo/dM)/{Jo/M) as defined bv iDemircan et al. I (|2006r i. Inserting 6 into (3) 



dP 5 dJo dM da 3d Jo , .dM 

-P=^^~6^1j; = - '^1^' T = (2 - ^)-^ = (26 - 3)—. (6) 

Change in mass ratio could be negligible because of a comparable relative mass loss from the 
component stars. Even ii dq/q 0, the term (1 — g)/(H-g) could be very small especially for 
high mass ratio (g ~ 1) systems. Thus, ignoring it in the first approximation is acceptable. 

Consequently, in eq. (6), 6 > 5/3 is required to decrease the period. But, 6 > 3/2 
is sufficient to shrink an orbit. If 3/2 < 6 < 5/3, orbital size decreases despite period is 
increasing. The size and the period of an orbit both increase i f 6 < 3/2. 



Using the mean decreasing rates of 0AM and mass from iDemircan et al. I 1)20061 ) , the 



mean value for the dynamical parameter (5) for detached CAB systems can be estimated as 



'77 _ Xdt _ -3-48 X 10^^° 

iM ~ dM ~ 1 on y 1 n-lO 
IT Mdt -^■'i^ X -LU 



S=^ = 1T^= ;on,,.n-in =2-68 (7) 



in the solar neighborhood. 

There may be alternative AM loss mechanisms among the detached CAB. Direct loss 
of 0AM is always possible in the galactic space due to stellar encounters, although it 
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would be negligible for short period systems (jStepien lll995HGliez. Neugebauer fc Matthews 



19931 ). Close encounters are more likely in multiple systems during periastron passages 



of distant third companions, which can lead to a rapid AM los s from the smaller orbits 



( Kiseleva. Eggleton fc Mikkila 



ferred to use 6 = A = 2.68 of 



1998 



Eggleton fc Kiseleva 1120011 1. Therefore, we have pre- 



Demircan et al. 



( 20061 ) as the mean dynamical parameter, 



which refers to all possible 0AM loss mechanisms when describing the mean dynamical 
evolution on log Jo — log P, log M — log P, log Jo — log M diagrams in the following sections. 



4 DISCUSSIONS 

4.1 The log Jo — logP and logM — logP diagrams 

The orbital angular momenta (Jo) and the periods (P) of CAB (Table 1) and W UMa (Table 
2) systems are all plotted on a log Jo — logP diagram in Fig. 1. CAB stars containing G, 
SG and MS designations and A & W types of W UMa stars are indicated. With smaller 
orbital periods (P <1 day) W UMa stars are concentrated at the lower left while CAB on 
the right display a wider band elongated towards the upper right. The constant total mass 
(M = Ml + M2) lines are computed using 



where P is varied while a chosen M is fixed. To represent a typical CAB, q = 0.88 median 
value for the present sample were used. 

The G systems with the biggest total masses prefer longer orbital periods and thus their 
OAMs are larger. The MS systems with smallest masses and shortest orbital periods have 
comparable OAMs with W UMa systems. Having moderate masses and orbital periods, the 
SG systems hold moderate OAMs. Unlike CAB with S, SG and MS designations, A and W 
type W UMa stars do not have distinct locations since they appear totally mixed. 

A well defined smooth upper boundary of CAB stars appears as if tracing a path of 
dynamical evolution. 0AM loss, mass loss and associated orbital period decrease would 
move a system form the upper right to the lower left parallel to the upper boundary. Finally 
some systems would enter in the r egion of con t act bi naries. With a similar idea and similar 



data set of less number of stars, 



Demircan 



( 19991 ) has computed dJo/dP directly from 



the inclination of the upper boundary line and estimated dP/dM and dJo/dM from its 
cut positions by the constant total mass lines. Decreasing of the total masses on the upper 
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boundary had been interpreted as the mass loss which activates 0AM loss. However, a logical 
reason why a system moves parallel to the upper boundary is not clear. The log Jo — log P 
diagram gives no clues on the time derivatives {dJo/dt, dP/dt and dM/dt) either. 

It is easy to explain why lower right of CAB distribution on Fig. 1 is empty by selection 
effects. Being less bright , small-mass, long-period systems could be missed by observers or 
have not yet studied. Such an excuse, however, does not exist to explain absence of short- 
period but more massive CAB systems on the diagram (upper left of CAB distribution). 
Therefore, we are inclined to think such systems do not exist at all. Otherwise, they would 
have been noticed and observed as some of the CAB stars in our list. 

In order to investigate why there exist such a well defined upper boundary, the log M — 
\ogP and log Jo — log P diagrams are compared in Fig. 2. Similar distribution characteristics 
are observable on both diagrams that the smooth upper boundary also exist on the log M — 
log P diagram. First, the upper boundary of CAB sample on the log M — log P diagram was 
eye estimated and digitized by computer. The estimated equation 

logM = 0.155 log P + 0.399 (9) 

seems to mark upper mass limits for various orbital periods. 

Next, from the P and M values of this line, we have computed the corresponding line 
on the log Jo — logP diagram using equation (8) with q = 0.88. Equally well, perhaps better 
fitting of the re-produced line on the log Jo — log P diagram clearly declares that the upper 
boundaries of these two diagrams are not independent and appears to be determined by mass 
upper limits for various orbital periods. Could it be that the mass upper limits imply initial 
Roche lobe structures with sizes primarily depend on the sizes of the preliminary orbits, 
where the mass within the lobes become limited according to the density of the star forming 
regions? or, is it only due to mass-period-activity relation? Why is the upper boundary not 
parallel to the constant total mass lines? (see Fig. 1). Special investigations other than this 
study seem to be needed. 

If period decrease occurs because of 0AM loss but no mass loss and transfer, a dynamical 
evolution would follow a path parallel to the constant total mass lines in Fig. 1, which would 
carry CAB stars into the empty region. Either this is not happening or systems moving into 
the region seize chromospheric activity. Therefore, the CAB upper boundary on log Jo — log P 
diagram may indicate a dynamical evolution with a minimu m mass loss. Neve r theles s, mean 



decreases determined from the kinematical ages of CAB bv iDemircan et al. I (|2006r i. which 
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are shown by the right sides of the triangles in both diagrams in Fig. 2, indicate that the 
direction of the mean dynamical evolution neither is parallel nor towards the empty region 
as it is shown by the arrows at the upper ends of the upper boundaries in both diagrams in 
Fig. 2. 

Both diagrams in Fig. 2 have been plotted intentionally to indicate MG and field sys- 
tems, in order see aging effect. Unlike the sub-groups with G, SG and MS designations, 
which prefers certain locations, kinematically young (ages < 0.6 Gyrs) MG systems appear 
randomly mixed among the kinematically old (3.86 Gyrs) field CAB. Random distribution 
of young and old systems in the CAB region confirms that an active binary can be born 
or may start its dynamical evolution anywhere on those diagrams. That is, a same location 
could be occupied equally likely by an old or a young system. Since there is no definite 
starting reference point (note: the direction and the speed may be known), it is not possible 
to trace the dynamical evolution of binary orbits on those plots in a similar manner as a 
single star nuclear evolution is traced on the H-R diagram. 

Unlike, MG and field CAB, which has no preferred location on the log Jo — logP and 
log M — log P diagrams, the W UMa stars designated with MG and field in Fig. 2 seem to 
show distinct locations which will be discussed in the following. 

4.2 logM — logP and log Jo — logP distributions of W UMa systems 

W UMa systems do not display similar distribution characteristics as CAB systems. Down- 
ward curvature of the upper boundary at the left end of Fig. 1 (W UMa region) imphes a 
decrease of 0AM because of hiding considerable fraction of total AM as a spin AM since 
synchronous rotation rates increase towards the shorter periods. Displaying the W UMa re- 
gion in a larger scale, Fig. 3 allows to compare mass and 0AM distributions. Not only 0AM 
distribution but also total mass distribution shows a similar trend. Therefore, the downward 
curvature of the upper boundary of W UMa region cannot be entirely due to transforming 
considerable 0AM to spin AM. It is more hkely to be related to mass content within the 
Roche lobes as the log M — log P diagram suggests. At smaller orbital periods, orbits are 
smaller. Smaller Roche lobes contain less mass. One should also not forget 0AM is more 
sensitive to orbital sizes than the total masses and periods as displayed in equation (1). 

The upper boundary of CAB region was extrapolated into the W UMa region and dis- 
played as solid and dashed lines in Fig. 3. One can easily notice that there are considerable 
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number of contact systems over those lines especially at upper right. Those are the systems 
with excess mass not only with respect to the rest of the W UMa sample but also with re- 
spect to the mass-period distribution inferred from the CAB systems. With an excess mass 
for a given orbital period, it is easier to fill Roche lobes in the early stages of formation. 

There are 18 systems above the solid line in Fig. 3a and half of them are with MG 
designations. Stars with MG designations have relatively larger masses and longer orbital 
periods with respect to the rest of the W UMa sample. Another interesting feature of W 
UMa distribution on Fig. 3a is that there is a square shaped empty region on the lower right 
corner that longer period systems (logP > —0.35) all have total masses bigger than 1.9 Mq 
(logM > 0.28). 

Displaying an unusual distribution with respect to the rest of the W UMa sample, the 
26 systems framed in a rectangle in Fig. 3 took our attention as young systems which were 
possibly born as contact binaries (with no pre-contact detached phases). We have restudied 
the galactic space velocities of those framed 26 systems, inspected their U-V diagrams, and 
computed their kinematical ages. Having the dispersions of 26.10, 19.16, 19.18 kms~^ at U, 



V, W and 37.63 kms^"*^ at the combined total space velocity vector S = \J{U'^ + V"^ + W^), 
where f/, V", W are galactic space velocity components with respect to LSR (Local Standard 
of Rest), an average age of 2.47 Gyr were found for them. It is also interesting that 15 out 
of 26 have positive V velocities bigger than Sun's V velocity in the LSR, which constitutes 
an additional strong argument favoring them to be young. Only two systems, e Cra and 
V2388 Oph, appear kinematically different. If these two systems were excluded, the average 
kinematical age reduces to 2.00 Gyr. 

Major contr ibution to the QAM loss at the contact stage comes from the gravitational 



wave radiation ( Guinan fc Bradstreet 



19881 ). which is inversely proportional to the orbital 



period. Having the largest orbital periods with respect to the rest of the W UMa sample, the 
systems in the frame are expected to be less effe cted by gravi t y wav es thus; their lifetime on 



the contact stage must be longer. According to 



Bilir et al. 



( 20051 ) ■ average lifetime of the 



contact stage is 1.61 Gyr, which is already s l ightly longer than the lifetime (0.1 < tcontact < 



1 Gyr) according to iGuinan &: Bradstreet I (|1988[ ). Thus, the 2.00 Gyr age found for the 
systems in the rectangle strongly favors them as contact binaries which were born contact. 
Existence of a flat lower boundary for the framed systems could be a hint for this. The 
systems below the flat boundary (systems with masses less than I.9M0), apparently do 
not had sufficient masses to be formed as contact systems. The systems below must have 
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been formed as detached systems and after their first Roche lobe overflow, their orbital 
periods decreased quickly due to mass transfer and they are now located among the other 
less massive contact systems. 

In fact, the region is not really empty but scarcely populated by detached systems (see 
Fig. 1 and Fig. 2). We believe those are the detached systems still evolving dynamically from 
the detached stage to the contact stage by loosing mass and orbital angular momentum. 
Before reaching a state of full contact, a mass transfer from more massive primary to less 
massive secondary must occur first. Such a mass transfer is not commonly observed because 
it is a fast stage at which the system mass and angular momentum appears to be conserved 
while the period of the orbit decreases quickly. If it is not possible to reach a full contact 
configuration in this stage, a detached binary may lose a great deal of its chance to be a 
contact binary and may become a classical Algol because mass loosing primary becomes 
less massive, then it is recognized as the mass loosing secondary. Mass transfer from the 
secondary to primary is a rather slow process which forces orbital period and size to be 
bigger. Only if there is sufficient 0AM loss, the orbital period may turn to decreasing. 

W type W UMa systems dominate over A types towards the shorter periods (see Fig. 
3). Therefore, according to the mass transfer evolution of an orbit as described above, W 
systems seems to be the first formed at the first mass transfer stage. Then, during the slow 
mass transfer stage after the mass ratio reversal, their orbital periods start to increase slowly 
and systems appear moving towards the longer periods into the region o f A type systems. 



Existence of many W systems with increasing orbital periods 



Oian 



( 20031 ) seems to c onfirm 



this scenario. A ty pes b eing formed after W types have been suggested by 



Mochnacki 



(19811), and 



Awadalla fc Hanna 



Wilson 



( 19781 ). 



( 200,51) . 



Just recently. 



Gazeas fc Niarchos 



( 20061 ). who studied mass-period distribution of over 



100 W UMa stars, have noticed A type systems usually have larger tot al masses. It is not 



Yakut fc Eggleton 



very o bvious but one can see Fig. 3b may also confirm this. According to 
( 20051). near contact bina r ies ca n transform to A type W UMa stars without a need of mass 

( 20061 ) propose A type systems may transform to W type systems 



loss. 



Gazeas fc Niarchos 



with a simultaneous mass l oss. The opposite were cla imed unreasonable since it requires 



the total masses to increase ( Gazeas fc Niarchos 



20061 ). However, it is clear in Fig. 3b that 



there are some A type systems with total masses less than the total masses of some W 
type systems. For a massive W system, there still could be chance to transform to a less 
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massive A type system. The two ways of forming A type systems cannot be ruled out only 
by comparing their mean masses. 

Do both A type and W type contact systems have a common origin? Is the one type the 
progenitor of the other? Such questions are still not answered. Formation of W UMa systems 
will be considered in the following section again when dynamical evolution from detached 
to contact state becomes clearer. 



4.3 The log Jo — logM diagram 

Because 0AM ( J^) and mass (M) are basic physical quantities determining orbital size and 
period (a & P), and because 0AM loss and mass loss are physical parameters controlling 
the magnitude and direction of dynamical evolution, the log Jo — log M diagram is a nat- 
ural choice to study dynamical evolution of binary orbits. Once, the diagram (Fig. 4) is 
produced, a sharp separation between the detached and contact systems stroke to our at- 
tention. Goodness of the separation is out striking that despite crowding along the border, 
there are only two systems (00 Aql, 6 Cap) on the wrong side, which could be due to a 
wrong identification of the state of being contact or just because of observational errors. 
Similar separation does not occur on the diagrams discussed before. 

Marking several positions on the borderline between CAB and W UMa stars, the follow- 
ing quadratic equation was produced. 

logJum = 0.522(logM)2 + 1.664(logM) + 51.315, (10) 

where M is in solar units and Jam is in cgs. Physical significance of this line is that it marks 
the maximum 0AM for a contact system to survive. It is like in single stars, spin AM has 
to be less than a certain value otherwise gravity cannot hold stellar mass together. If 0AM 
of a contact system is more than Jum, the contact configuration brakes. 

Because W UMa systems all have circular orbits (CAB systems of small orbital periods 
too) and because 0AM (eq. 1) is more sensitive to component separations than periods and 
masses, the line determined by (10) can be called "the contact border" . This is also because if 
it were possible to increase (Jo) of a contact system, there will be a limiting value associated 
with a limiting separation of components. Any Jo bigger than Jum, the distance between the 
components would be bigger than the limiting distance permitting a full contact. 

CAB svstems (RT Lac, AR Mon, e UMi, RV Lib, BH CVn) eliminated from the list of 
( 20061 ) ■ since they are filling or about to fill one of the Roche lobes that 
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mass transfer possibly occurring in them, are marked on Fig. 5. Since those systems are not 
close to the contact border and scattered randomly all over in the detached region, semi- 
contact configuration, like the other detached states, does not guarantee being close to the 
contact border. 

If total AM rather than 0AM were plotted on Fig. 4, the smoothness of the border 
would have been spoiled because as discussed by iRasio I (|l995l ) for low q contact binaries, 
the total spin AM of the primaries may be comparable with 0AM, so points of low q systems 
would be lifted by up to 1/3 of its value. It would have been meaningless to search such a 
border on a diagram using total AM. Moreover, adding total spin AM to 0AM would have 
introduced additional uncertainties since spin AM is less certain than 0AM because radii 
of gyration are uncertain. So it is not a coincidence for us to notice the contact border on 
the log Jo — log P diagram. 

Other diagrams (log Jo — log P, log M — log P) do not have such a well defined region and 
borderline for the contact binaries because according to equation (2) additional parameter 
(g) is needed to transform JumiM) function into Piim{M) function as 

Pum = 0.0046^^^t^Ml-5661ogM-0.008^ (11) 

where Pum is in units of days if M is in solar mass. It is also possible to express limiting 
separation of components as 

aum = 0.1167^^^^Mi-°^4i°s^^^+o-328, (12) 

where aum is in the units of solar radius if M is in solar mass. Because limiting periods (or 
separations) strongly depends on (g), a unique transformation of the contact border into the 
other diagrams is not possible. As displayed in Fig. 1 and Fig. 2, there are detached CAB 
systems with periods less than one day occupying the same region with W UMa stars on 
both log Jo — logP and logM — logP diagrams. With its own q, each system has its own 
Piim{M). Therefore, mixing CAB and W UMa stars is unavoidable on those diagrams. 

Being able to confine detached and contact systems into separate regions, log Jo — log M 
diagram is best for studying the dynamical evolution of binary orbits. The diagram clearly 
displays that there is a well defined possibility for a detached system to go into the region 
of contacts by losing 0AM. This possibility, however, depends on systems position on the 
diagram as well as speed and direction of its dynamical evolution, which could be very 
different from one system to another. This diagram too, like other diagrams, gives no clue 
on individual dynamical evolutions since there is no information about initial positions. 
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Random mixing of young (MG) and old (field) systems also occurs in this diagram. That is, 
a same position could be belong to both a young and an old system. 

Nevertheless, the contact border (eq. 10) could be helpful to determine contact binary 
candidates. The amount of mean losses (AJ and AM) corresponding to 2, 4 and 6 Gyr are 
subtracted from the Jo and M values of the border. That is, the contact border is shifted 
accordingly. Dotted lines in Fig. 5 represent shifted borders that the systems between a 
dotted line and the contact border are the ones, which have a definite chance to reach at the 
contact border within the time intervals indicated if their nuclear evolutions permit them 
to live as much. After counting, it becomes clear that 11, 23 and 39 per cent of the current 
sample of CAB could pass over the contact border within the next 2, 4 and 6 Gyr according 
to the mean dynamical evolution with 6 = 2.68. 

The dynamical evolution of contact binaries is more uncertain than detached CAB stars 
since a mean 6 is not available for them. If there is a smooth transition from the detached to 
the contact region, W systems are first to form according to W and A type distributions on 
the log Jo — log M diagram (see Fig. 4) because W systems are usually closer to the border 



than A type 



according to 



svstems. However, 



Gazeas fc Niarchos 



b rming first W types and then A types is not reasonable 
J ()2nn(Tl ) since this would need a mass gain. If A types are 
first to form, a smooth transition into the contact region also becomes unreasonable. This 
is because, forming W types later than A types appears also problematic since this require 
GAM gain (see Fig. 4). We encourage theoretical studies to investigate possibilities of jump- 
ing from the detached region into the contact region to form A type W UMa stars, which 
may occur during the first mass transfer stage, and then to produce W type W UMa stars. 
A to W or W to A transitions, if occuring, must be done within a lifet ime of a contact binary 



which is known to be very sh ort (0.1 < trnntart < 1 Gyr according to 
( 19881) . 1.61 Gyr according to iBilir et al. I ^W^)). 



Guinan Brad street 



4.4 Period and size evolution of orbits 

It is possible to draw constant period lines on the log Jo — logM diagram using (8). On the 
other hand, 

Jo = —^VGM^ (13) 

can be used to compute constant orbital size lines similarly. From the statistics of present 
CAB and W UMa samples, the median values of g = 0.88 and q = 0.39 are found and used 
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to represent CAB and W UMa systems in computing constant period and size lines which 
are shown in Fig. 6. The constant orbital period and orbital size lines run almost parallel 
to the contact border and both P & a values decrease towards it. Further decrease into the 
region of contacts is also clear. 

Constant period and constant size lines are sensitive to small q values. One can feel the 
sensitivity by comparing the lines of P = 1 day (or a = 3Rq) computed by g = 0.88 and 
q = 0.39, and P = 0.2 days (or a = 1.2Rq) computed with q = 0.39 and q = 0.1. Decreasing 
periods (or sizes) towards the lower right is deceptive as if evolution to contact stage is 
occurring from upper left to lower right which is impossible since such an evolution requires 
a mass gain. Mass loss, however, changes the direction from vertically down (0AM loss only) 
towards to the lower left (if 0AM and mass both are lost). Because mass loss dominates 
over 0AM loss (eq. 3), there are lower limits; one for orbital periods and one for the orbital 
sizes. Both limits are indicated by the dotted lines in the triangles and the corresponding 
numerical values (Alog Jo/AlogM) in Fig. 6. Any dynamical evolution with a 6 smaller 
than those limits indicates an increase rather than a decrease on both P a. 



5 CONCLUSIONS 

Total mass M and orbital angular momentum 0AM (Jo) are two basic physical quantities 
which determine sizes and periods of binary orbits. 0AM has several different analytical 
formulae established with different combinations of mass ratio g, semi-major axis a, orbital 
period P and systemic mass M. Because M h Jo are well known basic quantities and mass 
loss {dM) & 0AM loss {dJo) are again free parameters, from which the speed and direction 
of dynamical evolution can be inferred, the log Jo — log M diagram is the best for studying 
the dynamical evolution of binary orbits. 

There are limits in orbital dynamics, log Jo /(i log M > 3/2 is for decreasing orbital 
sizes, (i log Jo/ c? log M > 5/3 is for decreasing orbital periods, if d\ogJo/d\ogM takes a 
value between 3/2 and 5/3, orbital sizes shrink despite corresponding orbital periods in- 
crease. According to mea.n dyna mical evolution with 5 = dlogJo/dlogM = 2.68, as pro- 



vided bv lDemircan et al. I (|2006r i. the orbital sizes and periods of detached CAB systems are 
decreasing. 

A contact border, which separates W UMa stars from CAB, is discovered on the log Jo — 
logM diagram. So, One can predict 0AM of a system on the contact border, called Jum, 
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which is determined only by M. Any 0AM less than Jum implies a full contact for a binary. 
Similar to Jiim, the limiting period Pum, and limiting orbital size aum are computable using 
empirical formulae (11) and (12). Unlike Jum, which is independent of q, Pum and aum have 
a strong dependence on q. 

Rather than computing the time left for each detached system to reach at the contact 
border, limiting borders corresponding 2, 4, and 6 Gyrs were computed by subtracting mean 
decreases (A J and AM) of mean dynamical evolution from J and M values of the contact 
border. The CAB systems located between the contact border and those borders are the 
ones which are expected to be contact binaries within the time scales indicated. 

A detached system heading towards the contact region must go trough mass transfer 
stages before becoming a contact binary. At any mass transfer stage, even if it is conservative, 
0AM Jo is not conserved. For shrinking orbits, Jo must decrease according to eq. (1). For 
conservative cases, the 0AM difference between the two different size orbits of same system 
are to be compensated by spin AM. That is, Jo still decreases while the total spin AM is 
increasing. If a detach system is sufficiently close to the contact border, during the first fast 
mass transfer stage, it has a changes to jump over into the region of contacts and to be one 
of them. 

The traditional view of forming W UMa stars assumes that they only come from detached 
binaries of comparable periods. But, statistics available is convincing enough that there are 
very few such progenitors even to acc ount for the low s pace d ensity (0.2 pe r cent) of contact 



binaries in the solar neighbourhood (|Rucinski 



2002, 



20061) . Righteously 



Paczvnski et al. 



( 20061 ) says "at this time the contact systems seem to ap pear OUT OF NOWHE RE" in 

( 20061) . there- 



Paczvnski et al. 



order to emphasize the insufficiency of the traditional view, 
fore, consider mechanical tree-body orbital evolution, w hich is a similar formation model 
of close binaries in globular clusters (jPoolev et al. 1120031 ). where the inner binary has rela- 
tively longer orbital period and has a better chance to evolve into a contact system after all 

(20061) found that up to 50 per cent of W UMa binaries have compan- 



PribuUa fc Rucinski 



ions. 



The mean dynamical parameter 6 used in this study naturally contains affects of all 
mechanisms of 0AM loss since its value is derived from stellar kinematics. Perhaps, it is 
not possible now to make a firm conclusion which mechanis m is more efficient before a 



complete statistics and their analysis are available. Nevertheless, 



Bilir et al. 



20051 ) provided 



observational evidence that near contact detached binaries of few days orbital period are not 



Dynamical evolution of active detached binaries and contact binary formation 17 

the only source to form W UMa stars. They may also be born directly as contact systems. 
Moreover, it has also been found in this study that considerable fraction of current detached 
CAB sample (39 per cent within 6 Gyrs) to have a potential to jump over the contact 
border on the log Jo — log M diagram and become contact systems within a time scale of 
main-sequence life time. So, we encourage statistical studies to include these new findings 
after all contact binaries are rare with the local space density of just 0.2 per cent of the main 



sequence stars fjRucinski 



2002, 



200' 
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Figure 1. Orbital AM and period distributions of CAB and W UMa stars arc compared. Constant total mass lines (dashed) 
were computed using q=0.88, which is the median value for the present sample of CAB. 



Yakut K., Eggleton P.P., 2005, ApJ, 629, 1055 



20 Eker et al. 




V\/UIVla:l\/IG{ )) Field {•) 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

-1 -0.5 0.5 1 1.5 2 2.5 3 3.5 

Log P (days) 

Figure 2. Total mass-period (a) and OAM-Period (b) distributions are compared. CAB upper boundaries in both diagrams 
are not independent since it is possible to compute one from the other analytically. Young (MG) and old (field) CAB systems 
mix randomly. Mean dynamical evolution (arrow) and mean decreases (the right sides of triangles) are for 2 Gyrs. 
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Figure 3. W UMa regions of Fig. 1 and Fig. 2 are zoomed. Lines are CAB upper boundaries extended into W UMa region. 
Dashed is computed from solid using q = 0.39 (median value of W UMa sample). Young kinematical age (2 Gyr) and the empty 
region below imply systems in the rectangle may be formed directly as contact binaries (born in contact). 



22 Eker et al. 

54 



53.5 



53 



& 52.5 



^ 52 



51.5 



51 



50.5 



♦ ♦ + # ++ ^ 

.+ 




-0.4 -0.2 



0.2 0.4 

Log M(M®) 



0.6 0.8 1 



Detached region 



52.2 



52 



§) 51.8 



^ 51.6 



51.4 



51.2 




51 







Contact region 



0.1 0.2 0.3 0.4 0.5 
Log M(M®) 



Figure 4. A well defined borderline sharply separates deta<;hed and contact systems. Crowded region above (framed) is zoomed 
below. Symbols are like Fig. 1. 




Figure 5. Borders of equal times (dashed) to reach at the conta<;t border (solid). Mean dynamical evolution (hypotenuse) and 
mean decreases (right sides) for 2 Gyrs. 
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Figure 6. Constant period and size lines indicating Period and size evolutions of orbits. Mean (hypotenuse) evolution and 
corresponding decreases (right sides of triangles). Evolution with 5 = 5/3 (dashed) keeps P constant but evolution vi^ith (5 = 3/2 
(dashed) keeps a constant. 



Dynamical evolution of active detached binaries and contact binary formation 

Table 1. Physical pajrameters of the CABs. First 53 are possible members of young MG (< 0.6 Gyr) and the rest are 
(3.86 Gyr) field systems. 
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G5IV-III+A6V 


G 


3.95 





879 


71.614 


6.90 


1.70 


53 


101 





000 


47 


BD +17 703 


G2V+G8V 


MS 


1.99 





881 


75.683 


1.00 


1.00 


52 


614 





000 


48 


BM Cam 


KOIII 


G 


1.70 





545 


80.910 






62 


472 






49 


5 Get 


K2III 


G 


2.50 





786 


96.383 






52 


810 






50 


a Aur 


GlIII+KOIII 


G 


5.09 





951 


103.992 


12.80 


8.70 


53 


341 





000 


51 


V1817 Cyg 


K2III-II+A0V 


G 


7.73 





600 


108.893 






53 


623 






52 


T7 And 


G8III/IV 


G 


4.93 





904 


115.611 


11.00 


11.00 


53 


333 





000 


53 


SAO 23511 


F9.5V+G0V 


MS 


1.86 





682 


122.180 


0.75 


1.00 


52 


616 





000 


54 


XY UMa 


G0V+K5V 


MS 


1.748 





606 


0.479 


1.16 


0.63 


51 


762 





022 


55 


BI Cot 


G6V/IV+G6V/IV 


MS 


1.840 





916 


0.515 


0.90 


0.90 


51 


836 





024 


56 


SV Cam 


F5V+K0V 


MS 


1.786 





644 


0.593 


0.76 


1.18 


51 


815 





022 


57 


UV Leo 


GOV+G2V 


MS 


2.214 





970 


0.600 


1.19 


1.08 


51 


992 





019 


58 


BH Vir 


F8V-IV+G2V 


SG 


2.021 





981 


0.817 


1.11 


1.26 


51 


971 





014 


59 


WY Cnc 


G0-8V+K2? 


MS 


1.601 





495 


0.830 


0.58 


0.93 


51 


752 





014 


60 


CM Dra 


M4V+M4V 


MS 


0.444 





926 


1.268 


0.23 


0.26 


50 


938 





006 


61 


CC Eri 


K7.5V+M3.5V 


MS 


0.876 





537 


1.563 


0.41 


0.64 


51 


419 





Oil 


62 


V837 Tau 


G2V+K5V 


MS 


1.670 





673 


1.932 


0.74 


1.00 


51 


941 





004 


63 


EI Eri 


G5IV 


SG 


1.060 





379 


1.945 






51 


531 






64 


AR Lac 


G2IV+K0IV 


SG 


2.390 





897 


1.982 


2.72 


1.52 


52 


220 





008 


65 


BK Psc 


K5V+M4V 


MS 


1.040 





552 


2.168 


0.45 


0.72 


51 


594 





010 
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Table 1 - continued 



ID 


Name 


Spectrel 
Type 


Stage of 
Evolution 






1 


P 
(days) 


R\ 
(J) \ 


Jt2 


log Jo 

(cgs) 


J sp 


in! Jo 


66 


BH CVn 


F2IV+K2IV 


SG 


2.271 


0. 




2.612 


3.27 


3.10 


52.184 


0. 


,008 


67 


CF Tuc 


G0V+K4IV 


SG 


2.266 


0. 


.880 


2.799 


4.60 


1.50 


52.231 


0. 


.007 


68 


V711 Tau 


K1IV+G5IV 


SG 


2.531 


0. 


.821 


2.838 


3.80 


1.76 


62.310 


0. 


.008 


69 


PW Her 


K0IV+F8-G2 


SG 


2.652 


0. 


.768 


2.884 


3.80 


1.40 


62.343 


0. 


.007 


70 


AD Cap 


G5-8IV-V+G5 


SG 


1.624 


0, 


.526 


2.958 






61.966 






71 


TY Pyx 


G5V+G5-8V 


MS 


2.416 





.987 


3. 


199 


1.86 


1.58 


62.298 


0. 


.003 


72 


V1396 Cve 


M2V4-M4VC 


MS 


0.658 





696 


3.273 






51.346 






73 


HZ Aqr 


K3Vc+K7Vc 


MS 


1.234 


(] 


,804 


3. 


758 


0.45 


0.55 


51.830 





()()2 


74 


SZ Psc 


K1IV+F8IV 


SG 


2.861 





.766 


3.963 


5.10 


1.50 


52.444 


0. 


,006 


75 


Z Her 


K0IV+F5 


SG 


2.864 


0. 


.843 


3.990 


2.73 


1.85 


52.460 


0. 


.003 


76 


SAO 240653 


GOV+GOV 


MS 


2.180 


0. 


.998 


4.236 






62.264 






77 


UZ Lib 


K0III+A8? 


G 


1.440 





.309 


4.764 


21.00 


1.00 


61.839 





.019 


78 


RT Lac 


G5V+G9IV 


SG 


2.082 


0. 


.405 


5.070 


4.81 


4.41 


62.171 


0. 


.007 


79 


AS Dra 


G4V+G9V 


MS 


1.426 





.889 


5.408 






51.991 






80 


V1423 Aql 


G5V+G5V 


MS 


2.040 


1 


.000 


5. 


433 






52.252 






81 


BY Dra 


K6Vc + K7V 


MS 


0.934 


0, 


.891 




970 






51.699 






82 


RS UMi 


GOV+G-KV 


MS 


2.484 


0, 


.981 


6. 


166 






52.413 






83 


KT Peg 


G2V+K5V 


MS 


1.545 


0. 


.671 


6.209 


0.60 


1.00 


52.053 


0. 


,001 


84 


UX Ari 


A2/3V+K1/2V 


SG 


2.050 


0. 


.864 


6.442 


5.78 


1.11 


62.278 


0. 


.002 


85 


II Peg 


K2IV+M0/3V 


SG 


1.200 





.500 


6. 


730 






61.848 






86 


LR Hya 


K0/1V+K1/2V 


MS 


1.090 


0. 


.997 


6.871 


0.80 


0.80 


61.833 


0. 


.001 


87 


OU Gem 


K2V+K5V 


MS 


1.180 


0. 


.831 


6. 


998 






61.889 






88 


SS Boo 


GOV+KOIV 


SG 


1.928 


0, 


.988 


7.603 


3.30 


1.30 


52.260 


0, 


,001 


89 


MM Her 


G2IV+K1V 


SG 


2.469 





944 


7.962 


2.89 


1.56 


52.445 





OOl 


90 


FF Aqr 


G8III-IV+sdOB 


G 


3.100 





,240 


9.205 


6.00 


0.15 


52.427 





,0{)6 


91 


RIJ Cnc 


G8IV+Ffi-7 


SG 


2.930 





993 


10 


.163 


4.90 


1.90 


52.6(35 





,002 


92 


CQ Aur 


G8IV+F5V 


G 


3.764 





,882 


10 


.617 


9.91 


1.93 


52.791 


0. 


,003 


93 


RV Lib 


G8IV-K3IV 


SG 


2.785 


0. 


.183 


10 


.716 






62.294 






94 


EZ Peg 


G5IV+K0IV 


SG 


1.860 





.991 


11 


.668 






62.296 






95 


42 Cap 


G2IV+G2V 


SG 


2.375 


0. 


.727 


13 


.183 






62.480 






96 


AR Psc 


K1IV+G7V 


SG 


2.036 


0. 


.818 


14 


.289 


1.50 


1.50 


62.387 


0. 


,000 


97 


TZ Tri 


K0III+F5 


G 


5.099 


0, 


,976 


14 


.723 






53.060 






98 


V350 Lac 


K2IV~III 


G 


1.800 


0, 


,818 


17 


.742 






52.329 






99 


zcta And 


KlIII 


G 


3.480 


0, 


,289 


17 


.783 


13.40 


0.70 


52.653 





,006 


100 


UV CrB 


K2III 


G 


1.320 


0, 


,362 


18 


.664 






52.008 






101 


BL CVn 


KlII+FIV 


G 


2.672 





.991 


18 


.707 


15.20 


3.00 


52.627 


0. 


,003 


102 


AR Mon 


K2III+G8III 


G 


3.413 





.306 


21 


.232 


10.80 


14.20 


62.678 


0. 


.006 


103 


FG UMa 


G9III 


G 


2.080 





.387 


21 


.380 






62.370 






104 


IS Vir 


K2III 


G 


2.160 


0, 


.440 


23 


.669 






62.436 






105 


XX Tri 


KOIII 


G 


2.200 


0, 


.222 


23 


.988 






62.297 






106 


AI Phe 


F7V+K0IV 


SG 


2.440 


0, 


,968 


24 


.604 


2.93 


1.82 


52.601 


0, 


.000 


107 


IM Peg 


K2III-II 


G 


2.300 


0, 


,533 


24 


.660 






52.516 






108 


CS Get 


(G8-K1)III/IV+F 


G 


2.780 


0, 


,878 


27 


.353 






52.709 






109 


V792 Her 


F3V+K0III 


G 


2.881 


0, 


,960 


27 


.542 


12.80 


2.68 


52.737 


0, 


.002 


110 


TW Lep 


F6IV+K2III 


G 


1.990 


0, 


.951 


28 


.314 






52.473 






111 


V1762 Cyg 


K2IV-III+G8V 


G 


2.810 


0, 


.892 


28 


.676 


6.20 


0.90 


62.723 


0, 


.001 


112 


V965 Sco 


F2IV+K1III 


G 


3.418 


0. 


.990 


30, 


.974 


14.00 


6.60 


62.878 


0, 


.001 


113 


RZ Bri 


F0IV+G5-8III 


G 


3.323 


0, 


.963 


39 


.266 


6.94 


2.84 


62.892 


0, 


.000 


114 


V4200 Ser 


K2-3V+K2-3V 


MS 


1.700 


0, 


.998 


46 


.774 


0.80 


0.80 


62.432 


0, 


.000 


115 


EL Eri 


G8III-IV 


G 


1.930 


0, 


,379 


48 


.306 






52.430 






116 


AY Get 


WD+G5IIIO 


c; 


2.640 





,263 


56 


.885 


15.00 


0.012 


52.598 





,002 


117 


DK Dra 


KIIII+KIIII 


G 


3.480 


0, 


,981 


64 


.417 


14.00 


14.00 


52.997 


0, 


,001 


118 


V1197 Ori 


K4III 


G 


1.611 


0, 


,789 


142.889 






52.549 






119 


BD +44 801 


G2III+F2V 


G 


2.780 


0, 


.853 


963.829 






53.223 
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Table 2. Physical parameters of the W UMas. First 26 are possible members of young MG (< 0.6 Gyr) and the rest are older 
(5.47 Gyr) field systems. 



ID 


Star name 


Spectrel 


Type 


Mtot 




1 




p 


Ri 


R2 


log Jo 


Jap 


in/Jo 






Type 




(Mq) 






(days) 


(«0) 


(Rq) 


(cgs) 






1 


44 Boo 


K2V 


W 


1 


.339 


0, 


,487 





.2678 





.795 





.601 


61. 


,467 





.055 


2 


GZ And 


G5V 


W 


1 


.625 





,514 





.3050 


1 


.032 


1 


.032 


61. 


,624 





.056 


3 


TW Get 


G5V 


W 


1 


.964 





,530 





.3117 


1 


.048 





.819 


51 


,768 





.044 


4 


SW Lac 


KOVvar 


W 


1 


.738 





851 





.3207 





.997 





.925 


51 


724 





.047 


5 


RZ Gom 


F8V 


W 


1 


.592 


0, 


,437 


0, 


.3385 


1 


.088 





.728 


51 


,599 


0, 


.044 


6 


AG Boo 


F8Vn 


W 


1 


.968 


0, 


.403 


0, 


.3524 


1 


.314 


0, 


.572 


61. 


.744 


0, 


.046 


7 


V829 Her 


G2V 


W 


1 


.795 


0, 


.408 


0, 


.3582 


1 


.068 


0, 


.711 


61. 


.682 


0, 


.041 


8 


AH Gnc 


F5V 


W 


1 


.830 


0, 


,419 


0, 


.3604 


1 


.266 


0, 


.860 


61, 


,701 


0, 


.047 


9 


AE Phe 


F8V 


W 


1 


.995 


0, 


,461 


0, 


.3624 


1 


.269 


0, 


.843 


61, 


,781 


0, 


.042 


10 


AM Leo 


F5V 


w 


2, 


.009 





,449 





.3658 


1 


.260 





.983 


51 


,783 





.041 


11 


YY GrB 


F8V 


A 


1 


.777 


0, 


,243 





.3766 


1 


.426 





.805 


51 


,565 





.062 


12 


TX Gnc 


F8V 


w 


1 


.212 


0, 


,535 





.3829 


1 


.019 


0, 


.793 


51 


,450 


0, 


.037 


13 


BI GVn 


F2V 


A 


2, 


.325 


0, 


.413 


0, 


.3842 


1 


.200 


1 


.130 


61. 


.881 


0, 


.037 


14 


SS Ari 


GOV 


w 


1 


.749 


0, 


.302 


0, 


.4060 


1 


.361 


0, 


.795 


61, 


.618 


0, 


.047 


15 


AH Vir 


G8V 


W 


1 


.772 


0, 


.303 


0, 


.4075 


1 


.401 


0, 


.871 


61. 


.629 


0, 


.048 


16 


QX And 


F5V 


A 


1 


.412 


0, 


,261 


0, 


.4118 


1 


.632 


0, 


.804 


61. 


.430 


0, 


.078 


17 


UX Eri 


F9V 


A 


1 


,876 





373 





.4453 


1 


.362 





.678 


51 


728 





.038 


18 


NN Vir 


FO/FIV 


A 


2 


.(i7(i 





491 





.4807 


1 


.717 


1 


.246 


52 


()44 





.037 


19 


XZ Leo 


A8V 


A 


2 


.551 





,348 





.4877 


1 


.482 


1 


.283 


51 


950 





.038 


20 


DN Cam 


F2V 


W 


2 


.720 


0, 


,421 





.4983 


1 


.775 


1 


.224 


52 


()36 





.040 


21 


V351 Peg 


A8V 


W 


3, 


.165 


0, 


.360 


0, 


.5933 


1 


.875 


1 


.192 


52, 


,141 


0, 


.037 


22 


UZ Leo 


A9V 


A 


2. 


.703 


0. 


.303 


0, 


.6180 


2 


.024 


0, 


.851 


61. 


.996 


0, 


.043 


23 


BD +146016 


FOV 


A 


1 


.918 


0, 


.263 


0, 


.6369 


2 


.076 


1 


.181 


61, 


.707 


0, 


.052 


24 


V753 Men 


A8V 


W 


3, 


.296 


0, 


.970 


0, 


.6770 


1 


.738 


1 


.592 


62, 


.298 


0, 


.019 


25 


FN Gam 


A9V 


A 


2, 


.934 


0, 


.222 


0, 


.6771 










61. 


.988 






26 


II UMa 


F5III 


A 


2 


.623 


0, 


,172 





.8250 










51 


,861 






27 


CC Com 


K5V 


W 


1 


.050 


0, 


,522 





.2211 





.669 





.507 


51 


,263 





.068 


28 


V523 Cas 


K4V 


W 


0, 


.950 


0, 


,501 





.2337 





.695 





.522 


51 


,193 





.065 


29 


RW Gom 


KOV 


W 


1 


.230 


0, 


.337 


0, 


.2373 





.717 


0, 


.441 


51. 


,310 


0, 


.081 


30 


VZ Pbc 


K6V 


A 


1 


.510 


0, 


.911 


0, 


.2613 





.776 


0, 


.746 


61. 


.694 


0, 


.060 


31 


VW Gep 


G9V 


W 


1 


.144 


0. 


.276 


0, 


.2783 


1 


.209 


0, 


.398 


61. 


.234 


0, 


.097 


32 


BX Peg 


G4.5 


W 


1 


.400 


0, 


.373 


0, 


.2804 





.966 


0, 


.617 


61, 


.449 


0, 


.064 


33 


XY Leo 


KOV 


W 


1 


.305 


0, 


.600 


0, 


.2841 





.899 


0, 


.654 


61. 


.461 


0, 


.053 


34 


RW Dor 


KIV 


w 


1 


.070 


0, 


672 





.2854 





.823 





.636 


51, 


,342 





.041 


35 


BW Dra 


F8V 


w 


1 


.141 





281 





.2923 





.963 





.542 


51 


,244 





.083 


36 


OU Ser 


F9/G0V 


A 


1 


.194 





,173 





.2968 





.918 





.426 


51 


,145 





.113 


37 


TZ Boo 


G2V 


A 





.887 





,133 





.2972 


1 


.050 





.318 


50 


,846 





.167 


38 


FU Dra 


F8V 


w 


1 


.461 


0, 


,251 


0, 


.3067 


1 


.117 


0, 


.606 


51, 


,402 


0, 


.077 


39 


TY Boo 


G5V 


W 


1 


.330 


0. 


.437 


0, 


.3171 


1 


.072 


0, 


.761 


61. 


.469 


0. 


.054 


40 


YY Eri 


G5v 


w 


2, 


.160 


0, 


.403 


0, 


.3215 


1 


.163 


0, 


.793 


61, 


.798 


0, 


.045 


41 


FG Hya 


G2V 


A 


1 


.571 


0, 


.112 


0, 


.3278 


1 


.371 


0, 


.670 


61. 


.216 


0, 


.144 


42 


AO Gam 


GOV 


w 


1 


.680 


0, 


.413 


0, 


.3299 


1 


.042 


0, 


.817 


61, 


.624 


0, 


.043 


43 


AB And 


G8V 


W 


1 


.499 


0, 


,560 





.3319 


1 


.041 





.755 


51 


.589 





.042 


44 


W UMa 


F8V 


w 


1 


.760 


0, 


,479 





.3336 


1 


.126 





.813 


51 


,684 





.048 


45 


EQ Tau 


G2V 


A 


1 


.754 


0, 


,442 





.3413 


1 


.139 





.786 


51 


,672 





.048 


46 


VW Boo 


G5V 


W 


1 


.400 


0, 


,428 





.3422 










51 


,504 






47 


V757 Gen 


F9V 


W 


1 


.690 


0, 


.690 


0, 


.3432 


1 


.054 


0, 


.880 


51, 


,701 


0, 


.044 


48 


V508 Oph 


F9V 


A 


1 


.530 


0. 


.616 


0, 


.3448 


1 


.049 


0, 


.811 


61. 


.698 


0, 


.048 


49 


V781 Tau 


GOV 


W 


1 


.738 


0, 


.406 


0, 


.3449 


1 


.166 


0, 


.776 


61. 


.661 


0. 


.042 


50 


BT Leo 


G8V 


w 


1 


.120 


0, 


.342 


0, 


.3465 










61, 


.301 






51 


BV Dra 


F7V 


W 


1 


.398 


0, 


.402 


0, 


.3501 


1 


.101 


0, 


.742 


61, 


.496 


0, 


.050 


52 


GK Boo 


F8/7V 


A 


1 


.986 


0, 


111 





.3552 


1 


.097 





.815 


61, 


.394 


0, 


.096 


53 


QW Gem 


F8V 


W 


1 


.700 





,334 





.3581 


1 


.242 





.924 


51, 


,602 





.050 


54 


BB Peg 


F5V 


W 


1 


.946 





,360 





.3615 


1 


.286 





.833 


51 


,717 





.046 


55 


LS Del 


GOV 


W 


1 


.467 


0, 


,375 





.3638 


1 


.046 





.800 


51 


,522 





.047 


56 


V410 Aur 


G0/2V 


A 


1 


.412 


0, 


,144 


0, 


.3663 


1 


.397 


0, 


.605 


51, 


,239 


0, 


.121 


57 


V752 Gen 


F7V 


W 


1 


.730 


0, 


.311 


0, 


.3702 


1 


.280 


0, 


.754 


61. 


.604 


0, 


.049 


58 


V417 Aql 


G2V 


W 


1 


.900 


0, 


.362 


0, 


.3703 


1 


.290 


0, 


.825 


61, 


.706 


0, 


.045 


59 


XY Boo 


F5V 


A 


1 


.081 


0. 


.168 


0, 


.3706 


1 


.246 


0, 


.566 


61. 


.077 


0. 


.113 


60 


HV Aqr 


F5V 


A 


1 


.564 


0, 


.146 


0, 


.3745 










61, 


.319 






61 


U Peg 


G2V 


W 


1 


.528 


0, 


,330 





.3748 


1 


.207 


0, 


.735 


61, 


.629 


0, 


.051 


62 


RT Lmi 


F7V 


A 


1 


.774 


0, 


,367 





.3749 


1 


.272 


0, 


.838 


61, 


.669 


0, 


.043 


63 


HX UMa 


F4V 


A 


1 


.720 


0, 


,291 





.3792 










51, 


,688 






64 


EE Get 


F2V 


W 


1 


.829 


0, 


,315 





.3799 


1 


.350 





.820 


51 


,651 





.046 


65 


AU Ser 


KOV 


A 


1 


.567 


0, 


,701 


0, 


.3865 


1 


.051 


0, 


.950 


51, 


,665 


0, 


.039 
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Table 2 - continued 



ID 


Star name 


Spectrel 
Type 


Type 


Mtot 




1 


P 
(days) 


Rl 
(«©) 


H2 
(«0) 


log Jo 

(cgs) 


J sp 


in/Jc 


66 


BH Cas 


K4V 


W 


1. 


.096 


0. 


,474 


0. 


.4069 


1 


.114 


0, 


.751 


61 


.367 


0, 


.037 


67 


HT Vir 


F8V 


A 


2, 


.280 


0. 


812 


0, 


.4077 


1 


.262 


1 


.252 


61 


.963 


0, 


.034 


68 


EX Leo 


F6V 


A 


1, 


.866 


0. 


199 


0, 


.4086 










61 


.667 






69 


V839 Oph 


F7V 


A 


2, 


,136 


0. 


305 


0, 


,4090 


1 


.462 


0, 


.989 


61 


.766 


0, 


.047 


70 


V566 Oph 


F4V 


A 


1 


826 


0. 


243 


0, 


,4096 


1 


.515 





.769 


51 


.597 





.060 


71 


UV Lyn 


F6V 


W 


1 


844 


0, 


,367 





4150 


1 


.387 





.889 


51 


.702 





.040 


72 


RZ Tau 


A7V 


A 


2 


,516 


0. 


540 





,4157 


1 


.594 


1 


.055 


51 


.992 


0, 


.040 


73 


BL Eri 


GOV 


W 


0, 


,939 


0. 


542 





,4169 


1 


.021 


0, 


.777 


51 


.279 


0, 


.031 


74 


V842 Her 


F9V 


W 


1, 


.713 


0. 


260 


0, 


,4190 










51 


.571 






76 


Y Sex 


F8 


A 


1, 


.430 


0. 


,182 


0. 


.4198 


1 


.240 


0, 


.620 


61 


.341 


0, 


.073 


76 


V899 Her 


FBV 


A 


2, 


.876 


0. 


566 


0, 


.4212 


1 


.504 


1 


.171 


62 


.096 


0, 


.036 


77 


AK Her 


FBV 


A 


1, 


.610 


0. 


229 


0, 


.4216 


1 


.657 


0, 


.820 


61 


.494 


0, 


.077 


78 


ER Ori 


F7V 


W 


2, 


.160 


0. 


552 


0, 


.4234 


1 


.388 


1 


.130 


61 


.883 


0, 


.035 


79 


EF Dra 


F9V 


A 


2 


,102 


0. 


160 


0, 


,4240 


1 


.704 





.779 


51 


.582 





.097 


80 


EF Boo 


F5V 


W 


2 


,349 


0. 


512 


0, 


,4295 


1 


.090 


1 


.460 


51 


.939 





.027 


81 


AP Leo 


F8V 


A 


1 


,916 


0. 


297 


0, 


,4304 


1 


.507 





.838 


51 


.689 





.047 


82 


AW UMa 


F2V 


A 


1 


,370 


0. 


070 





,4387 


1 


.692 





.560 


50 


.988 





.233 


83 


V776 Cas 


F2V 


A 


1. 


,948 


0. 


130 


0, 


,4404 


1 


.710 


0, 


.710 


51 


.465 


0, 


.113 


84 


TV Mus 


F8V 


A 


1. 


.473 


0. 


119 


0, 


.4467 


1 


.643 


0, 


.742 


61 


.236 


0, 


.130 


86 


VB02 Oph 


GOV 


W 


1, 


.778 


0. 


,335 


0. 


.4634 


1 


.493 


0, 


.935 


61 


.670 


0, 


.043 


86 


AA UMa 


F9V 


W 


2, 


.192 


0. 


545 


0, 


.4680 


1 


.356 


1 


.101 


61 


.910 


0, 


.029 


87 


DK Cyg 


A8V 


A 


2, 


.343 


0. 


,325 


0, 


.4707 


1 


.789 


0, 


.987 


61 


.868 


0, 


.051 


88 


V728 Her 


F3V 


W 


1 


,949 


0, 


178 





4713 


1 


.789 





.898 


51 


.576 





.082 


89 


EL Aqr 


F3V 


A 


1 


,880 


0, 


203 





4814 


1 


.730 





.880 


51 


.591 





.069 


90 


AH Aur 


F7V 


A 


1 


,967 


0, 


169 





.4943 


1 


.856 





.897 


51 


.573 





.085 


91 


DO Aql 


G2V 


A 


1 


,920 


0, 


,846 


0, 


.5068 


1 


.382 


1 


.283 


51 


.862 





.029 


92 


V401 Cyg 


FOV 


A 


2, 


,166 


0. 


,290 


0, 


,5827 


1 


.950 


1 


.170 


51 


.816 


0, 


.046 


93 


e Cra 


F2V 


A 


1. 


.940 


0. 


,128 


0. 


.6914 


2, 


.200 


0, 


.788 


61 


.600 


0, 


.109 


94 


AQ Tuc 


F2/5 


A 


2, 


.620 


0. 


,358 


0, 


.6948 


2, 


.027 


1 


.296 


62 


.003 


0, 


.041 


96 


RR Cen 


A9V 


A 


2, 


.243 


0. 


,210 


0. 


.6067 


2, 


.188 


0, 


.974 


61 


.762 


0, 


.066 


96 


V535 Ara 


A8V 


A 


1, 


.980 


0. 


,303 


0, 


.6293 


1 


.847 


1 


.121 


61 


.772 


0, 


.038 


97 


AG Vir 


A8V 


A 


2 


,120 


0, 


317 


0, 


,6427 










51 


.835 






98 


S Ant 


A9V 


A 


2 


700 


0, 


,392 


0, 


,6483 


2, 


.070 


1 


.360 


52 


.056 





.035 


99 


V1073 Cyg 


F2V 


A 


2 


,110 


0, 


,319 


0, 


,7859 


2 


.154 


1 


.318 


51 


.862 





.034 


100 


V2388 Oph 


F3V 


A 


1 


,954 


0, 


186 


0, 


,8023 


2 


.373 


1 


.368 


51 


.668 





.060 


101 


TY Pup 


F3V 


A 


2, 


,920 


0. 


183 


0, 


,8192 


1 


.787 


1 


.082 


61 


.957 


0, 


.035 


102 


V921 Her 


A7IV 


A 


1, 


.974 


0. 


226 


0, 


.8774 










61 


.744 







